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o 

^ i The non-linear effects operating at the recombination epoch generate a non-Gaussian signal in the 
, CMB anisotropics. Such a contribution is relevant because it represents a major part of the second- 
order radiation transfer function which must be determined in order to have a complete control 
■ of both the primordial and non-primordial part of non-Gaussianity in the CMB anisotropics. We 
' provide an estimate of the level of non-Gaussianity in the CMB arising from the recombination epoch 
OO . which shows up mainly in the equilateral configuration. We find that it causes a contamination to 
^) . the possible measurement of the equilateral primordial bispectrum shifting the minimum detectable 

value of the non-Gaussian parameter Z^*^'' by A/^*j^'' = 0(10) for an experiment like Planck. 
^ ■ 
OO 

o 

1 Introduction 

b 

I Cosmological inflation pj has become the dominant paradigm to understand the initial conditions for 
the Cosmic Microwave Background (CMB) anisotropics and structure formation. This picture has 
recently received further spectacular confirmation by the Wilkinson Microwave Anisotropy Probe 
(WMAP) five year set of data [2]. Present [2] and future [3] experiments may be sensitive to the 
non-linearities of the cosmological perturbations at the level of second- or higher-order perturbation 
theory. The detection of these non-linearities through the non-Gaussianity (NG) in the CMB Q has 
become one of the primary experimental targets. 

A possible source of NG could be primordial in origin, being specific to a particular mechanism 
for the generation of the cosmological perturbations. This is what makes a positive detection of 
NG so relevant: it might help in discriminating among competing scenarios which otherwise might 
be undistinguishable. Indeed, various models of inflation, firmly rooted in modern particle physics 
theory, predict a significant amount of primordial NG generated either during or immediately after 
inflation when the comoving curvature perturbation becomes constant on super-horizon scales [4]. 
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While single-field [5J and two(multi)-field [6] models of inflation generically predict a tiny level of 
NG, 'curvaton-type models', in which a significant contribution to the curvature perturbation is 
generated after the end of slow-roll inflation by the perturbation in a field which has a negligible 
effect on inflation, may predict a high level of NG [TJ [8]. Alternatives to the curvaton model are 
those models characterized by the curvature perturbation being generated by an inhomogeneity in 
the decay rate [HI [TO] or the mass [llj or of the particles responsible for the reheating after inflation. 
Other opportunities for generating the curvature perturbation occur at the end of inflation [12] and 
during preheating [13]. All these models generate a level of NG which is local as the NG part of 
the primordial curvature perturbation is a local function of the Gaussian part, being generated on 
superhorizon scales. In momentum space, the three point function, or bispectrum, arising from 
the local NG is dominated by the so-called "squeezed" configuration, where one of the momenta is 
much smaller than the other two and it is parametrized by the non-linearity parameter /^l- Other 
models, such as DBI inflation |14j and ghost inflation [15], predict a different kind of primordial NG, 
called "equilateral", because the three-point function for this kind of NG is peaked on equilateral 
configurations, in which the lenghts of the three wavevectors forming a triangle in Fourier space are 
equal [TO]. The equilateral NG is parametrized by an amplitude /^l^'' [I3|- Present limits on NG 
are summarized by -9 < f^l < 111 and -151 < f^f^ < 253 at 95% CL jSUH]. 

On the other hand there exist many sources of NG in the CMB anisotropies beyond the primor- 
dial ones, which are essential to characterize in order to distinguish them from a possible primordial 
signal. One should account for the so-called secondary anisotropies, which arise after the last scat- 
tering epoch. For example, cross-correlations SZ-lensing and ISW-lensing [191 ISO] produce a bias 
in the estimate of NG which is at the level of the expected estimator variance at Planck angular 
resolution [21]. Analogous conclusions have been reached in Ref. [22j for the cross correlations of 
density and lensing magnification of radio and SZ point sources with the ISW effect. Furthermore, 
the impact of cosmological parameters' uncertainties on estimates of the primordial NG parameter 
in local and equilateral models of NG has been recently studied in [23] . 

There exists another relevant source of NG: the non-linear effects operating at the recombina- 
tion epoch. The dynamics at recombination is quite involved because all the non-linearities in the 
evolution of the baryon-photon fluid at recombination and the ones coming from general relativity 
should be accounted for. The first steps in describing the physics at recombination at second-order 
in perturbation theory were taken in |241 125j (see also Ref. [27\ I26j). where the full system of Boltz- 
mann equations at second-order describing the evolution of the photon, baryon and cold dark matter 
fluids were obtained (see also Ref. [28]) These equations allow to follow the time evolution of the 
CMB anisotropies at second-order at all angular scales from the early epoch, when the cosmological 
perturbations were generated, to the present through the recombination era. Such a contribution is 
so relevant because it represents a major part of the second-order radiation transfer function which 
must be determined in order to have a complete control of both the primordial and non-primordial 
part of NG in the CMB anisotropies and to gain from the theoretical side the same level of precision 
that could be reached experimentally in the near future [4j. 

The NG generated at the surface of last scattering comprises various effects, as described in 
details in Ref. [25] (for some specific effects see also [271 129]). It turns out that the dominant 
contribution comes from the non-linear evolution of the second-order gravitational potential which 
grows in time on small scales. Since this effect is a causal one, developing on small scales, we expect 
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that the NG it generates will be of the equilateral type, rather than of local type. As we will see our 
results confirm such an expectation. Therefore, a reasonable question is to which extent the NG from 
recombination alters the possibile detection of the primordial NG of the equilateral type. The goal of 
this paper is to estimate in a semi-analytical way the contribution to NG from recombination. Along 
the same lines of Ref. 130], we adopt a simple analytical model to parametrize the transfer functions 
and test its goodness in Section 2 by evaluating the minimum value of equilateral NG detectable 
by the Planck experiment for which there exist numerical calculations [31]. As a by-product of our 
results, we find that the signal-to-noise ratio for a primordial equilateral bispectrum scales as the 
square root of the maximum multipole £max probed by an experiment, unlike the well known scaling 
as ^max for the local case. In Section 3 we compute the bispectrum of the NG generated by the 
evolution of the second-order gravitational potentials, and in Section 4 we provide an estimate of 
the NG from recombination showing that it corresponds to a degradation in the measurement of an 
equilateral primordial bispectrum of A/^'l''' = 0(10), shifting the minimum detectable value from 
f^^^ — 67 to /^'l''' — 79 for an experiment like Planck. 

2 Signal-to-Noise ratio for the primordial equilateral 
bispectrum 



In this Section we wish to recover the estimate for the signal-to-noise ratio (S/N) given in Ref. |31j 
for the primordial bispectra of "equilateral" type [T7] by adopting a simple model. In other words, 
we test the goodness of the semi-analytical model we will be using in the next Section to estimate 
the bispectrum from the recombination era. 

Our starting point is the primordial equilateral bispectrum [T7] 

(<I>(ki)cI>(k2)<I>(k3)) = (2^)3^(3) (1^^ ^ ^ k3)i?equil(fcl, A;2, ^3) , (1) 

where 

pequil fi42'l 1 1 21 



5equii(fci, A:2, ks) = f^l • 6A • ( --3-3 - -3-3 - -3-3 - -2-2-2 + + (5 perm.) , (2) 

y /v-j^ (1^2 1 3 2 3 1 2 3 ixij^ nji2 iv^ j 

and the permutations act only on the last term in parentheses. The parameter /^'l''' quantifies the 
level of NG while A = 17.46 x 10~^ is the amplitude of the primordial gravitional potential power 
spectrum computed at first-order 

(cI>W(ki)$«(k2)) = (27r)V3)(ki + k2)P(fci), (3) 

with P{)i) = A/k"^ . Since the signal-to-noise ratio {S/N) will be some function of the maximum 
multipole a given experiment can reach, ^max ^ 1, we can use the flat-sky approximation j3Ul I32| 
and write for the bispectrum 

{a{ii)a{i2)a{i^)) = (27r)2<5(2)(423)i?(^i,^2,4) , (4) 
where fi23 = 4 + ^2 + ^3, with [T6] 



5equii(4,4,^3) = / ^^A;fdA:fdA:|5(i)(A:f23)i?cquii(M,^2,fc3)A^(4>^DA^(4,^l)A^(4,^": 



(2vr)2 



(5) 
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where k' means k evaluated such that fc" = ^/(tq — t^) and 



^0 dr 



A^{£, k^) = / ^5(V(fc-)2 + ^7(To-r)2, r)e^^'i-r-r) ^ ^g) 

Jo (^0 - t)^ 

is the radiation transfer function defined by the CMB source function S{k,T). In this notation, tq 
and Tr represent the present-day and the recombination conformal time, respectively and k^ and A;" 
are the momentum components perpendicular and parallel respectively to the plane orthogonal to 
the line-of-sight. 

The {S/N) ratio in the flat-sky formalism is [30^ I32j 

where /sky stands for the portion of the observed sky. In order to compute the bispectrum -Bequii(^i> ^2, ^3) 
and the power spectrum C{i) we adopt the following model 



— e^'='("°-"'-)$(k')A^(^,A;") (8) 
where we mimic the effects of the transfer function on small scales as 

A^{e, F) = a (ro - r,)-2e-i/2(£/^.)^-%-i/2(|fc.|A.)^-^ , (9) 

i.e. a simple exponential and a normalization coefficent a to be determined to match the amplitude 
of the angular power spectrum at the characteristic scale £ ~ = A;*(ro — Tr)0 It is important 
to make clear what are the reasons underlying the choice of such a model. When computing the 
(S/N), Eq. ([8]) with i^: = fc*(ro — Tr) ^ 750 and a ~ 3 is able to account for the combined effects 
of "radiation driving" , which occours at ^ > £eq — 160 and boosts the angular power spectrum with 
respect to the Sachs- Wolfe plateau, and the effects of Silk damping which tend to suppress the CMB 
anisotropies for scales i > in — 1300. The combination of these effects produces a decrease in the 
angular power spectrum from a scale £^ ~ 75oEI The power spectrum in the flat-sky approximation 
is given by {a{h)a{l2)) = {27rf6^^\h2)C{ei) with 

C{i)=^^^^^ I dk^\A^{£,kn\'P{k). (10) 

The exponential of the transfer function for Eq. ([8]) allows to cut off the integral for /c ~ A;* and one 
finds (see also Ref. [[30] ) 

where the last equality holds for i ^ i^. To compute the bispectrum we proceed in a similar 
way. One first uses the Dirac deltas, ^^^H^m) and d^'^^lus). Then it proves to be useful the 
change of variable kf = xiIi/{to — Tr), fcf = X212/{tq — Tr). In this way the transfer functions 



^We could equally choose a transfer function as A^{l,k^) = a{TQ — Tr) ^^^^ _ j/c^l)^ the 

relevant approximation being that the integral over k^ is cut at the scale k^. 

^ The choice of the exponent 1.2 derives from the study of the diffusion damping envelope in Ref. [35]. 



4 



become A'^{£i,k^) oc e"^/^^'^*'^'/^*^^'^ which allows to cut the integrals over Xi {i = 1,2) at 
Now, as a good approximation to see the effects of the transfer functions, we can take £ ^ and 
thus the integral over Xi can be easily computed by just evaluating the integrand in Xj = times 
4:{£^/ii){£^/£2)- With this approximation the integral in kf is easily obtained and we get for the 
bispectrum 

B.,^,{£^,£2,£s) = ^/e.-V^2^-(4-^+4'^+4-^)/2^P^2 



p3f3 f3p3 03 03 02p2n2 ' P^o2f3 



+ + (5 perm.) , (12 



where 



'"^ - "l + £l + 2£i-£2. (13) 



'^3 



The coefficient /i ~ 1/1.4 = 0.7 is a fudge factor that improves the matching between our ap- 
proximation for the bispectrum and numerical results that have been consistenly checked. Notice 
that, according to our approximation, the equilateral structure of Eq. ([2]) is preserved in £ spacejfl 
In computing the signal-to-noise ratio, consistency with our approximation ()12p requires that we 
integrate over £i , £2 starting from a minimum ^min > I* up to £max and paying attention to the fact 
that even ^3 in Eq. (fT3]) must be larger than ^min- The scaling with £yaa.x with respect to the case of 
a local type bispectrum turns out to be much milder. While for the local type {S/N)"^ oc ^^ax fSO] - 
for the equilateral bispectrum ([l]) we fincfl (S/N)"^ cx ^max and, setting £^ = 750 and ^min — 1200, 

I) ' = 0.48 X 10^ ^ ^ (/Z")' W - 8 /sk, A if^T'f ^max. (14) 
By choosing /g^y = 0.8 and ^max = 2000 we find a minimum detectable 

/Nf'-ee, (15) 

obtained imposing (5'/A^)equii = 1- Both the estimate of the minimum value of f^^^ and the 
scaling (S/N)'^ oc £max are in remarkable agreement with the result obtained in Ref. [31] where 
the full transfer function is used and a value of /^l^'' = 67 is obtained!^ Notice that our estimate 
is independent from the coefficient a and the exponetial e~^^'^^^^^*^^'^ introduced below Eq. ([8]) to 
mimic the full transfer function. This is because there is an equal number of transfer functions in 
the numerator and denominator of the expression ([Tj) for the signal-to-noise ratio and their effect 
tend to cancel despite they are not simple multiplicative factors (see discussion in Ref. 



■^The expression (jl2p can be also written as 
5oquii(^i,^2,^3)^= i^nrHSh /^q-ia3 A2£2e-(^}-^+4-^+4-^)/2^i-^ (1 + cos0)(4 +i2~ £3)/£l£¥l being the 
angle between £1 and £2- 

* These scalings can be easily understood by analyzing the expressions (S/N)'^ for the local and equilateral 
primordial NG. In the local case, (5/iV)2 is proportional to / (f£i(f£2d'^£3S'^^H£i23){el+il+il)^/{£i£2e3f [30]; 
since the squeezed configuration, e.g. £\ £2,(3, is dominating the local bispecrum, the integral becomes 
proportional to / d£id£2{£2 / £1) oc ^^^^x- t^e equilateral case, however, (S/N)^ receives contributions from 
the configuration which is peaked at £1 £2 ^ £3 and therefore it can be written as / (P£i(P£2S{£i — ^2)/^? oc 

-^-max- 

^Wc thank M. Liguori for discussions about the minimum value of /nl"' detectable by Planck and for its 
scaling with €max- 
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3 Non-Gaussianity from recombination 



Comforted by the goodness of our model, in this Section we wish to estimate the level of NG 
generated at the recombination era. One can check that on small scales the second-order anisotropics 
are dominated by the second-order gravitational potential ^^"^^ which grows as r^, as first pointed 
out in Ref. [25] Therefore, the main contribution to the bispectrum generated at recombination 
comes from 

eri = ^Ag + ^^^ (16) 

which is the usual term appearing in the CMB anisotropics due to the intrinsic photon energy density 
fluctuations Aqq and the gravitational redshift due to the potential [25] . Such a term on large scales 
reduces to the Sachs- Wolfe effect, while on small scales at recombination 

^Ag +cl>(2) -i?$(2) = -^G(ki,k2,k)r(A:i)cI>«(ki)T(A;2)<I>(i)(k2)r2, (17) 

where we have evaluated the expression at the recombination time Tr and R = Sp^/Ap^ is the 
baryon-to-photon energy density ratio. Eq. (17) is the extension to the second-order of a well-known 
expression at linear order [33] (see also [3ll EH]), as discussed in details in Ref. [26]. It can also be 
easily obtained using the expressions of Ref. [25] by taking the sound speed as a function of R. The 
kernel is given by [2S] 

G(k,,k.,k)=k,.k.-^iiiJ^l]|iJ^. (18) 



Notice that this expression has been obtained assuming all the momenta much larger than keq |25j . 
From the form of the kernel we see that the NG at recombination is dominated by an equilateral 
configuration, as expected from the fact that its origin is gravitational. Here and in the following 
we are implicitly assuming that a convolution is acting on the kernel as 

J dkidk2 <5(3)(ki + k2 + k3)G(ki,k2, k)r(fci)$(i)(ki)T(A;2)$(^)(k2) . (19) 



(2vr)- 

The reader should remember that, at first-order in perturbation theory, the combination Gg^+i? 
is exponentially suppressed by the Silk damping, but still greater than the term R^^^^ (which does 
not suffer the damping) for the maximum multipole of interest, £max ~ 2000. This is meanly due to 
the fact that the first-order gravitational potential rapidly decays on small scales. On the contrary, 
at second-order in perturbation theory, the gravitational potential grows like the scale factor on 
small scales and it turns out that the R^^^^ dominates on small scales (see Ref. [26j). 

The gravitational potential at linear order can be expressed as usual in terms of the transfer 
function T{k) 

T{k)^ 12 (^^^ ln[fc/8fceq], (20) 

where the last step is an approximation valid on scales smaller the the equivalence scale, k ^ keq. In 
the following we will account for the logarithmic growth just with a coeffcient ro(fc) = 12 ln[fc/8fceq] ~ 
11 for the scales of interest. 



^ Notice that we have numerically verified that such a term starts to be relevant from multipoles £ ~ 1300. 
This is actually indicated also by the results of Ref. [5^] . 
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In the flat-sky approximation one arrives at an expression similar to ([5]), where now one of the 
linear transfer functions must replaced by a transfer function at second-order. Specifically one finds 



rccl^l,«2,«3j 



(tq - TrY' 

(2vr)2 



dA;f(ifc|dA:|(^(i)(A:f23)[G(k;,k^,k'3)T(A;;)r(A:^)P(A:i)P(A:^) (21) 



X {ii,kf)A^ {£2, kl)A^'^^\e3,kl)+ cyclic 
By using our model ([8]) and 



A^^^\£,kn 



R 



14 (ro - Tr 



(22) 



for the second-order radiation transfer function, we find 

R (to - r,)-4 



reclti,t2,t3j 



A;,\r2A2a2r2e-i/2(£,/£,)^-^e-i/2(^2/^*)- / dkfdkldklS^'Hk 



123) 



1 



G'(ki,k2,k3) ^ e 



-l/2(|fci,|/fc.)i-2 _i/2(|fc2.|/A..,)l 



'^1 '^2 



+ cyclic 



(23) 



At this point we proceed further by employing the same approximation described after Eq. (jlip . We 
use the Dirac delta to replace the variable k^z, and the exponential allow us to evaluate the integral 
for kiz = k2z = 0, for scales £i » This leads to 



Srcc(4,^2,4) 



4/2 R ,2 



(27r)2 14 -^Ov'^cqTrj <^eq^*^ ^ 

10 (4-^"'i)(4-/2)- 



h-l2 



1^2 '- 



^2 

*-3 



+ cyclic . 



(24) 



Again here /2 is a coefficient to better calibrate our approximations with numerical results that we 
have performed in order to test the validity of our approach. Not surprisingly, it turns out that 
/2 ~ /i ~ 1/1.4. 



4 Contamination to primordial non-Gaussianity from 
recombination: Fisher matrices 

Our goal now is to quantify the level of NG coming from the recombination era and to estimate the 
level of degradation it causes on the possible measurement of the equilateral primordial bispectrum. 
The reader should keep in mind that, given the form of the kernel function (jlSp . the NG from 
recombination is expected to be of the equilateral type. A rigorous procedure is to define the Fisher 
matrix (see, for example, [20] ) 

where i (or j)= (rec,equil), and to define the signal-to-noise ratio for a component i, {S/N)i = 
\l \j F^^, and the degradation parameter di = FnF^^ due to the correlation bewteen the different 

components rij = Fl[^^ / J Fl[l^ F~-^ . The first entry -Feqml.equil of the Fisher matrix corresponds to the 
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{S/N)'^ ratio computed in Eq. (jl4p which does not account for any kind of cross-correlation. Due 
to the equilateral form of the NG generated at recombination we expect that the minimum value 
detectable for /^'l''' will be higher that the one reported in Eq. (jlSp . For the mixed entry we find 



rec,equil 



/sky 1 



vr (27r)2 

/sky 4i? 48 r2 



d«l(X«2d«3 toj 6C(^i)C(^2)C(£3) 



-3/1/2 



7r3 14 2^6 a 



1 



(^1+^2-^3) 



4-^2 



10 (^3-4)(^3-^2) 



(26) 



where £3 is given by Eq. (jl3p . The factor 3 in front of this expression comes from cyclic permutations. 
The integral can be performed numerically and, integrating from a minimum ijahi — 1200 up to 
•^max = 2000, and by taking R ~ 0.3 when evaluated at recombination, a ~ 3, Tq ~ 11, (/i:cq7"r)^ — 26, 



i-eq 



150, 



750, we find -^rec,equii — ^ 10 ^- Finally for the entry i*rec,rec W6 get 



TT (27r)2 



6C{h)C{e2)C{£3) 



= /: 

+ 6 



2 /sky f ^\ f^O \ ^ \4 d 

2 257r36 \ 14 / \ a 



cq 



1.2 £3 



1^2 



10(V_W3_ 

3 £i 



id£2ei/2(£3/^.)^-^ei/2(^2/^.) 

10 (/2-/3)(^l-^"'2)^ 



^f^2^i 



10(£3-^l)(^3-^2) 



^3 



^2 



and we find a value i^rec.rec — 0.014. We are now able to compute the entries of inverse of the Fisher 
matrix, -^jj ^- In the following we report our results for the signal-to-noise ratios and the degradation 
parameters 



^^/cquil 
S_ 

N 



'"rec,equil 



equil,equil 



12.6 X 10"^/, 



-3 ^equil 
NL ' 



0.1, 



F, 



rec,rec 
-1 



rcc,equil 



-0.53, 



equil 



equil, equil i'ec,rec 

F F^^ ~ 1 4 

reCirec-* roc, roc — ^-^ ) 

-fcquil, equil -f'cquil, equil — 1-4 • 



(28) 
(29) 

(30) 

(31) 
(32) 



As a confirmation of our expectations, we find that the NG of the type given by Eq. (jl8p has a quite 
high correlation with an equilateral primordial bispectrum. This translates into a degradation (or a 
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contamination) in the mimimum detectable value for f^'^^ with respect to the value given in (jlSp . 
In fact from the signal-to-noise ratio (j28p we find a minimum value of 

/Z" ^ 79 , (33) 

imposing that (57A^)cquii = 1- This roughly corresponds to a contamination to the primordial 
equilateral NG ofl 

A/X' = 0(10) . (34) 

A similar way to quantify this statement is by defining an "effective" f^l^ for which the equilateral 
bispectrum ([2]) has the same Fisher matrix errors as the recombination bispectrum (see also |36i [3T] ) 



NL 



free 
Jl 



/■cquil_ 



(35) 



\/-^equil,cquil 

In this case we find an effective non-linearity parameter f^f^ ~ 8 which agrees with the result (|34p l^ 
Similarly we can compute the Fisher matrix accounting for the NG generated at recombination and 
the primordial NG of the local type 

(<I>(ki)$(k2)$(k3)) = (27r)35(3) (ki + k2 + ka) Aoc(A:i, ^2, h) , (36) 

where 

BUh, h, ks) = f^l • • (lia + i + • (37) 

The bispectrum and the signal-to-noise ratio as defined in Eq. d?]) have already been computed in the 
flat-sky approximation in Ref. [30j. The result is that iS/N)f^^ = 47r-2/sky(^*/Cin)(/NL)^^Cax> 
corresponding to a minimum detectable value of = ^(7) for ^max = 2000 (when other possible 
sources of NG are ignored) . We can compute the off-diagonal entry of the Fisher matrix in a similar 
way to what we have described in this section, and we get -Frec,ioc — 8 x IO'^/^IJ'l- Finally the entry 



^Due to a non- vanishing correlation, ry, {S/N) gets modified from its zero-order value to (S/N) = 
(S'/iV)o(l-r2 )l/^ so that the minimum detectable value of /^l"'' gsts shifted by a quantity A/^^"''/ {f^'^^^)o = 

^ An alternative quantity can be used in order to measure the contamination to the primordial bispectra. 
It is that effective equilateral which minimizes the defined as 



One finds 



(/nL ^cq(^l, ^2, k)\f^^^^^i — -Brcc(^l, ^2, ^3, 
6C(4)C(^2)C(^3) 

-^rcCjCquil 



2 



/rec 

/nl — p 

-^^cquil.cquil 



and an analogous expression to compute the contamination to the local primordial bispectrum. In both cases 
we find a similar value to the one obtained from Eq. ([35]) and Eq. (|43)l . Notice however that the effective 
non-linearity parameter defined in this way contains a somewhat richer information with respect to (j35|) : we 
are not just comparing signal-to-noise ratios, but we are asking what is the value of equilateral (local) /nl 
which best mimics the bispectrum from recombination. 
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0.014 has already been computed above. From inverting the Fisher matrix, we get the 
following signal-to-noise ratios and the degradation parameters 

^ 14x10-2/;-, (38) 



s 



loc, loc 

:~0.1, (39) 



roc, loc n, A A f Ai^\ 

rrecioc = , _^ _ = ^ -0.44 , (40) 

Y -^loc.loc-^fec.rec 

dree — -^rec,rec-^rec,rec — 1-2, (41) 

C^loc = -^equil,equil-^equil,equil — ' (4^) 

In particular, from Eq. (j38p we see that now the minimum detectable value of remains basically 

unchanged in the presence of the recombination signal. Similarly the effective reads 



= Vgec^ -0.7, (43) 

Y -f^loc,loc JNL ^ 

which is much smaller than the effective non-linearity parameter (j35p for the equilateral case. We 
have also checked the cross-correlation between the primordial local and equilateral bispectra finding 
a value of rioc,equii — 0.23, which is in agreement with the value reported in Ref. [31]. This reflects 
the fact that the primordial local and equilateral signals are not fully uncorrelated. The reason is 
due to the fact that the equilateral and local bispectrum (j37p and ([2]) approach the same shape in the 
equilateral configuration. This is also the reason why the cross-correlation between the primordial 
local and recombination bispectra is not so small. 

We conclude with two comments. First, we would like to stress that the NG from recombination 
and the one due to the non-linear evolution of gravity from the last scattering surface to us include 
many other contributions. However, they will affect the measurement of the primordial local NG 
|37j . while the contamination to the primordial equilateral NG is dominated by the non-linearities 
considered in this note. Second, our estimates have been obtained with a maximum multipole 
of ^max = 2000. For an experiment like Planck, it approximately corresponds to the angular scales 
where the instrument noise and the secondary effects from lensing are still negligble when computing 
the signal-to- noise ratios [20l [30] . 
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